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It is shown that the idea that scaling behavior in turbu- 
lence is limited by one outer length L and one inner length 
r\ is untenable. Every n'th order correlation function of ve- 
locity differences F n [Ri, R2, • • ■) exhibits its own cross-over 
length r/ n to dissipative behavior as a function of, say, Ri. 
This length depends on n and on the remaining separations 

Ra , R3 , One result of this Letter is that when all these 

separations are of the same order R this length scales like 
r) n (R) ~ v(R/L) Xn with x n = (Cn - Cn+i + Cs - <a)/(2 - Ca), 
with £ n being the scaling exponent of the n'th order struc- 
ture function. We derive a class of scaling relations including 
the "bridge relation" for the scaling exponent of dissipation 
fluctuations /i = 2 — Qq. 
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The aim of this Letter is to expose the fact that the 
notion of the dissipative length in hydrodynamic turbu- 
lence is a rich and interesting concept whose complex- 
ity exceeds the expectations of established models and 
standard theories [JO. Indeed, during a few decades the 
thinking about the universal small scale structure of tur- 
bulence was dominated by Kolmogorov's picture of en- 
ergy cascade through an "inertial interval" which is lim- 
ited on one side by the integral scale of turbulence L 
and on the other side by the Kolmogorov viscous scale 
r\ = (i^/e) 1 / 4 where v and e are the fluid's kinematic 
viscosity and the mean energy flux in the turbulent flow 
respectively. During the last decade there has been a 
growing concern about the inability of Kolmogorov's the- 
ory to cope with the increasing experimental evidence for 
multiscaling (or multifractal) behaviour of higher order 
structure functions. Together with the concern about 
the statistical theory there arose a realization that the 
uniqueness of the viscous length is suspicious. Paladin 
and Vulpiani M , and also Frisch and Vergassola j| used 
the multifractal model of turbulence to assess the charac- 
teristic viscous lengths associated with the higher order 
structure functions of velocity differences 



S 2n (Ri) = <Kn|ri)H , 
s 2n+ i(Ri) = Ri • (uKnkDKnK)! 2 ") , 



(i) 

(2) 



where wfrxlr^t) = u(r[,t) — u(ri,t) and u(r,t) is the 
velocity field of the fluid, R\ = r\— r\, and R\ = R\/R\. 
In homogeneous and locally istoropic turbulence S n (R) is 
a function of the magnitude of R, and the viscous length 
is that value of R at which the functional dependence of 
S n (R) changes from a non-trivial power law S n (R) ~ R^ n 
to a trivial power law that stems from a Taylor expansion 



of the velocity differences, S n (R) ~ R n . The multifractal 
model leads to a prediction that this length depend on 
the order n. In this Letter we argue that a proper dis- 
cussion of cross-overs to dissipative behaviour requires 
the analysis of functions richer than structure functions. 
Firstly, we state that the fundamental object to analyze 
is the n-point correlation of velocity differences 

F n (r- L ,r' 1 ...r n ,r' n ) = (w(ri\r[) . . . w(r n \r' n )) , (3) 

which is an n-rank tensor. All the separations Ri = \r[ — 
Tj I and Tij = |r, — rj\ are within the "inertial range". It 
is generally accepted that this correlation function is a 
homogeneous function of its arguments, i.e. 

^„(An, Xr[ . . . Xr' n ) = X u T n {r x , r\ ...r' n ). (4) 

It should be understood that quantities like (Q) are ob- 
tained from (|^) by fusing some coordinates together. (In 
this case all — > and all Ri — > R). In this process of 
fusion one crosses the viscous scale, and it is important 
to understand how to do this. 

Our discssion will not call for any ad-hoc model of tur- 
bulence. It will be based on two solid building blocks, 
one being the Navier-Stokes equations, and the other the 
fusion rules that were derived recently. The fusion rules 
appear naturally in the analytic theory of Navier-Stokes 
turbulence j|-|7j and passive scalar turbulent advection 
, and they determine the analytic structure of the 
n-order correlation functions (^) when a group of coor- 
dinates tend towards each other. The fusion rules were 
derived in |?j for systems in which Eq. (|^) holds with uni- 
versal scaling exponents (i.e the scaling exponents do not 
depend on the detailed form of the driving of the tur- 
bulent flows). The fusion rules address the asymptotic 
properties of T n when a group (or groups) of coordi- 
nates tend towards a common coordinate withing each 
group, while all the other coordinates remain separated 
by a large distance R. There are two particular examples 
of fusion rules that we will employ in this Letter. The 
first pertains to the fusion of one pair of points. When 
the distance between one pair is small, R\ ~ p, and the 
separations between all the other coordinates are much 
larger, Ri ~ R for i ^ 1, then to leading order in p/R 



T n cx S n (R)S 2 (p)/S 2 (R). 



(5) 



The second situation pertains to the case in which we 
have two groups of fusing coordinates separated by a 
large distance R. When there is a group of p points 
separated by a typical distance p\ , and a group of n — p 
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points separated by a typical distance pi with a large 
distance R between the groups, then 

T n ~ S n {R)S p {pi)S n -p(p 2 )/Sp{R)S n -p(R). (6) 

These forms hold as long as p, p\ and pi are in the inertial 
range. 

The Navier-Stokes equations for an incompressible ve- 
locity field u(r, t) may be written in the form 

u(r, t) + P [u(r, t) • V] u(r, t) = vV 2 u(r, t) . (7) 

Here v is the kinematic viscosity and P is the transverse 
projector. Given the equation of motion we can take the 
time derivative of Eq. (||) . We find 

n 

jF n = £ (wirr \r[, t) . . . w(r, \r' p t) . . . w(r n \r' n ,t)). 

Substituting Eq. (0) , and considering the stationary state 
in which dlF n /dt = we find the balance equations 

X> n (n,ri; . ..r n ,r' n ) = J n (r x ,r' x ; . ..r n ,r' n ) . (8) 

The term J „ originates from the viscosity term in (]?]), 

n 

J n ( ri y i; . . . r n y n ) = uJ2 (V 2 + V 2 ,) (lotnlri) . . . 

...w(r j \r' j )...w(r n \r' n )) . (9) 

The term T> n stems from the nonlinear term, and it needs 
a bit of algebra to bring to the exact form 

/n 
dr^P^pir) (10) 

xiw^irxlr'J . ..L ( \r h r'^r) . . . w an (r n \r' n )) , 

n 

L?{r o y v r) = -^{^( r .- r | rfc )V7 (11) 

n k=i 

+w 1 (r' j -r\r' k )V],}w fi (r j - r\r'j - r). 

We are going to argue now that when all the separations 
Rj are of the same order of magnitude R, the interaction 
term has a very simple evaluation, i.e. 

V n ~S n+1 (R)/R. (12) 

To this aim we need to prove that the integral is local in 
the sense that it converges in the ultraviolet and in the 
infrared. 

As the coordinate r is being integrated over, the most 
dangerous ultraviolet contribution comes from the re- 
gion of small r. In this region the projection operator 
can be evaluated as 1/r 3 . Other coalescence events of 
r with other coordinates contribute less divergent inte- 
grands since the projection operator is not becoming sin- 
gular. When r becomes small, there are two possiblilites: 



(i) rj ^ r,t and (ii) rj — tv In the first case the corre- 
lation function itself is analytic in the region r — > 0, and 

we can expand it in a Taylor series Const + B ■ r + 

where B is an r-independent vector. The constant term 
is annihilated by the projection operator. The term lin- 
ear in r vanishes under the dr integration due to r — > — r 
symmetry. The next term which is proportional to r 2 is 
convergent in the ultraviolet. In the second case we have 
a velocity difference across the length r. Accordingly we 
need to use the fusion rule (||), and we learn that the 
leading contribution is proportional to r^ 2 . This is suf- 
ficient for convergence in the ultraviolet. We note that 
the derivative with respect to rj cannot be evaluated as 
1/r when rj — r^. Rather, it is evaluated as the inverse 
of the distance between rj and the nearest coordinate in 
the correlation function. 

FIG. 1. Typical geometry with (n — 1) velocity differences 
in a ball of radiusi? on the left separated by a large distance 
r S> R from a pair of points on the right. 

To understand the convergence of T) n when the in- 
tegration variable r becomes very large we consider the 
relevant geometry as shown in Fig.l. There is one veloc- 
ity difference across the coordinates rj — r and r'j — r 
(which is shown on the right of the figure), (n — 1) ve- 
locity differences across coordinates that are all within a 
ball of radius R (at the left of the figure) , and one veloc- 
ity difference across the large distance r which is much 
larger than R. In the notation of this figure the leading 
order contribution for large r is obtained from the fusion 
rules (|J) for the situation on the right and ([|) for the 
geometry on the left. The resulting evaluation for the 
leading term is r^ n+1 (Rj /r)^ 2 (R/r)^- 1 . On the face of 
it, this term is near dangerous. For any anomalous scal- 
ing the integral converges since Cn+i < Cn-i + C2 due to 
Hoelder inequalities. This convergence seems slow. How- 
ever, the situation is in fact much safer. If we take into 
account the precise form of the second-order structure 
function in the fusion rules we find that the divergence 
with respect to rj translates in fact to dS^ 1 (Rj) / dRj~ t 
which is zero due to incompressibility. The next order 
term is convergent even for simple (K41) scaling. This 
completes the proof of locality of ([l(]) . The conclusion is 
that the main contribution to the integral in ( [To| ) comes 
from the region r ~ R. Therefore the integral can be 
evaluated by straightforward power counting leading to 
(|i"2|). It should be stressed that a more detailed analy- 
sis demonstrates that when the separations between the 
coordinates that do not involve velocity differences, (i.e 
separations like rjk but not Rj)go to zero, the evaluation 
does not change. 

The evaluation of the quantity J '„ is more straight- 
forward. When all the separations Rj and r,-j are of the 
same order R, the correlator in (^|) is evaluated simply 
as S n (R). The Laplacian is then of the order of 1/R 2 . 
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We note that when v — > (which is the limit of infi- 
nite Reynolds number Re), this term becomes negligi- 
ble compared to D n . The ratio J n /T> n is evaluated as 
i/S n (R)/RS n +i(R), which for fixed R vanishes in the 
limit v — > 0. Thus the "balance equation" becomes a 
homogeneous integro-diffcrcntial equation T> n = which 
may have scale-invariant solutions with anomalous scal- 
ing exponents Cn+i ( n + It should be stressed 
that the evaluation jl^ ) remains correct for every term 
in 2?„, but various terms cancel to give zero in the homo- 
geneous equation, provided that the scaling exponent C, n 
is chosen correctly. To make this important point clear 
we exemplify it with the simple case n — 2 for which 
D n can be greatly simplified. Consider the scalar object 
J : 2{rx,T l 1 ,r 2 ,r' 2 ) = (w(ri\r[) ■ w(r 2 \r' 2 )). The terms in 
the scalar balance equation for this case are exactly 

V 2 (r 1 ,r[,r 2 ,r' 2 ) = d[S 3 (r 12 >) - S 3 {r 12 )]/2dr x 

+ d[S 3 (rv2) - Ss(ri'3')]/2dri , (13) 

J2(ri,r' 1 ,r 2 ,r' 2 ) = i/{V?[5 2 (ri 2 /) - S 2 (r 12 )} 

+ Vl[S 2 (ry 2 )-S 2 (r V2 ,)}} , (14) 

When all the separations are of the order of R we can see 
explicitly that J 2 ~ uS 2 (R) / ' R 2 which is much smaller 
than each term in D 2 . Considering the scale invariant 
solution S 3 (R) = AR& where A is a dimensional coeffi- 
cient, we see that 



12' 



r l'2 



Cs-X 
"l'2' 



Obviously the solution for T> 2 = requires the unique 
choice £3 = 1 which is the known exponent for S3 |Q. 
The coefficient A is now determined as e which is the 
mean energy dissipation per unit mass and unit time. 

There is a cross-over from the scale invariant solu- 
tion of the homogeneous equation to dissipative solutions 
when J 2 becomes comparable to any of the terms in T> 2 . 
This happens when at least one of the separations ap- 
pearing in ( pT| ) becomes small enough. Denoting the 
smallest separation as r m we evaluate J 2 ~ vS 2 (r m )/r 2 a . 
From this we can estimate, using the balance equation, 
S 2 (r m ) ~ (S 3 (R)/vR)r 2 a ~ er^jv. In the inertial range 
we have S 2 (r) ~ (er) 2 / 3 (r/L) ?2 ~ 2 / 3 . The viscous scale 
r] 2 for the second-order structure function is then deter- 
mined from finding where these two expressions are of 
the same order of magnitude, i.e. 



e V 2 2 /u=(e V2 ) 2 / 3 (r/L)^ 3 . 



Using the outer velocity scale Ul we estimate e ~ U^/L 
and end up with 



V2 



LRe-V(2-C 2 ) . 



(16) 



Note that this result is not in agreement with the ad-hoc 
application of the multifractal model jl]-[§ which predicts 
V2 ~LRe- 2 /( 2+ ^. 



A similar mechanism operates in the general case of 
n 2. As long as all the separations are in the iner- 
tial interval J~ n is negligible. When one separation e.g. 
ri2 diminishes towards zero, and all the other separations 
are of the order of R, the internal cancellations leading to 
the homogenous equation T>„ = disappear, and T> n is 
evaluated as in (|l2[). The term J7* n is now dominated by 
one contribution that can be written in short-hand nota- 
tion as v\I\F n (r 12 , {R}). We can solve for F n (ri 2 , {R}) 
in this limit: 



F n {r 12 ,{R}) 



r 2 2 S n+1 {R)/vR . 



(17) 



On the other hand we have, from the fusion rule (|^), 
the form of the same quantity when r\ 2 is still in the 
inertial range, i.e. F n (r 12 , {R}) sa S 2 (r 12 )S n (R)/ S 2 (R). 
To estimate the viscous scale r\ n we find when these two 
evaluations are of the same order. The answer is 



/ R 

r) n (R) = m{j j 



Cn + C3 — Cn+1 - C2 ,-. „, 
X n = ^ J • (18) 



We note that the Hoelder inequalities guarantee that 
x n > and increases with n. We see that the viscous 
"length" is actually an anomalous scaling function. 

Next we show that in the same spirit we can derive 
important (and exact) scaling relations between the ex- 
ponents Cn °f the structure functions and exponents in- 
volving correlations of the dissipation field. We consider 
correlations of the type 

Kt ] = (e(r 1 )w{r 2 \r' 2 )...w{r n \r' n ))^R-^\ (19) 
Kg = (e(rMr[)w(r 2 \r' 2 )...w(rnK))^R~^\ (20) 

where R is a typical separation between any pair and 
e(r) = i/|V«(r)| 2 , and we are interested in the scaling 

relations between the exponents fi n and the exponents 

(2) 

£ n - Note that /!q in this notation is the well studied 
Jl0| , p| exponent of dissipation fluctuation which is de- 
noted as fi. This relation is almost at hand for jin \ 
We see this by writing 



K^=v lim ViV 2 ^„(r 12 ,{i?}). 

ri2->0 



(21) 



Using the result ([l?]) we find immediately 

/iW=l-Cn+3- (22) 
(2) 

The scaling relations satisfied by ii n require considera- 
(15) tions of the second time derivative of the correlation (|^). 



K= (w(r 1 \r[,t)...w(r l \r' l ,t) 

j,i=x 

...w(r j \r' J ,t)...w(r n \r' n ,t)). 



(23) 



Using the Navier-Stokes equations for the time deriva- 
tives we derive a new balance equation T>^ = J {2) 
where, using the definition (Jllf), 
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/n 
drdr' £ P(r)P(r')(w(r 1 \r' 1 ) 



(24) 



'■.;=i 



.L(n,ri,r) . ..Lfa^'pr') . ..w(r n \r' } 



Using the fusion rules and following steps similar to those 
described above, we can prove that the integrals over r 
and r' converge. Accordingly, when all the separations 

(2) 

are of the order of i?, every term in Dh, ' is evaluated as 
S n+2 {R)/R 2 . The term J {2) takes on the form 



3 ( ?=v 2 £ (v? + v?,)(v?+v?0 



(25) 



x (w(ri\r[)...w(r i \r' i )...w(r j \r' j )...w(r n \r' n )) . 

As before, when all the separation in this quantity are 
of the order of i?, the Laplacian operators introduce 
factor of 1/R 2 and the evaluation of this quantity is 
J {2) ~ v 2 S n {R)/R A . Clearly this is negligible compared 



(2) ! 

to typical terms in D\ ■ The quantity 13^' contains a 
cross contribution with one Laplacian operator and one 
nonlinear term with a projection operator. The integral 
is again local, and one can show that the evaluation is 
~ vS n+ i(R)/ R 3 which is also negligible compared 

(2) 

to typical terms in D„ ■ 

Now we consider the fusion of two pairs of coordinate, 
e.g. J"i2 — > and r^ — > 0. As before, the cancellations in 
T>^ are eliminated, and the evaluation of a typical term 
becomes the evaluation of the quantity. The other two 
terms in the balance equation also become of the same 
order because the Laplacian operators Vf and V3 are 
evaluated as r^ 2 and r% 2 respectively. As before we can 
consider the resulting balance equation as a differential 
equation for F n {r\ 2 , ^34, {R}). The leading term in this 
equation is 

Av 2 V 2 V 2 F n {r 12 ,r^,{R}) » ~ S n+2 (R)/R 2 . 



(2) 



The solution is 



F n (ri2, r 3i , {R}) ~ r 2 12 r 2 3i S n+2 (R)/v 2 R 2 



(26) 



Finally we can write the quantities /Q™' in terms of the 
correlation function as 

K.£>=v 2 lim V 1 V 2 V 3 V 4 ^„(r 12 ,r 3 4,{i?}). (27) 

ri2,r-34— »0 

Using (pff) here we end up with the evaluation 



~S n+6 /R 2 ccR-^ 2 \ 



(2) _ 



2-Cn 



+(>■ 



(28) 



(2) 

For the standard exponent /1 = p/ a we choose n = 
and obtain the phenomcnologically proposed "bridge re- 
lation" fi = 2 — ( 6 . To our best knowledege this is the first 
solid derivation of this scaling relation. In general, if we 



have p dissipation fields correlated with n velocity differ- 
ences the scaling exponent can be found by considering 
p time derivatives of (||), with the final result 



M, ( t P) =P-(n+3p- 



(29) 



We see that Eqs.(|22j), (28) and (29) can be guessed if we 
assert that for the sake of scaling purposes the dissipa- 
tion field e(r) can be swapped in the correlation function 
with w 3 (ri|r' 1 )/i?i, where R\ is the characteristic scale. 
This reminds one of the Kolmogorov refined similarity 
hypothesis. We should stress that (i) our result does not 
depend on any uncontrolled hypothesis, and (ii) it does 
not imply the correctness of the hypothesis. Our result 
is implied by the refined similarity hypothesis, but not 
vice versa. 
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